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1. INTRODUCTION 

The electric supply industry has been changing rapidly throughout the world as a consequence of the 
implementation of deregulation. In deregulation, the user of the interconnection network has to pay the cost 
of transmission service, which is determine of transmission line usage each user [1]. Many methods have 
been used or proposed to evaluate of the transmission usage allocation (TUA) through tracing method. 

Bialek [2] proposed the principle of Proportional Sharing Method (PSM) to determine the 
contribution of each generator to the load. However, this method only considers positive contribution to the 
line flow. Meanwhile, Rudnick [3] proposed the method of tracking the power flow algorithm along with the 
Proportional Sharing principle by using a distribution method through the DC load flow model. Kirchen [4] 
proposed a tracking method using the principle of proportional assumption, which is a common method. 
This principle is based on the concept of grouping the buses and branches in the network system into 
homogeneous group. This method can also be used to calculate the losses contribution of a generator in the 
line. The disadvantages this method is that it cannot determine reverse power flow. F.Wu [5] developed a 
graph theory which can be used to calculate the contributions of each generator and loads to line flows, and 
also the real power transfer between individual generators and loads, while Yu [6] proposed a Tee 
Proportional Method (TPM) by using the principle of proportional development approach without using the 
inverse matrix. However, the method has not been tested for tracing reactive power. 

Gubina [7] proposed a new method for determining the generator's contribution to a particular load, 
which uses the nodal generation distribution factors (NGDF-s). Pantos [8] recommended the modifications of 
the factors in determining GGDF to prevent expansion of the matrix, by using the load distribution factor 
(LDF) and the load distribution factor (LSF). This method can also be used to overcome the capacitive 
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loading on the line. Valerie.S [9], on the other hand, proposed the concept of reference frame of the loop by 
modelling the network through injection of voltage and current at each node. All flow in the network is a 
representation of power delivery in a series of related loops. Hussain Shareef [10] suggested using modified 
nodal Equations to identify the real and reactive power transfer between generators and load. 

Meanwhile, S.Abdulkader [11] recommended a tracking method based on the flow with the actual 
physical line, where the contribution of the division of the load current is in accordance with the cross- 
sectional area of conductors. This concept is expressed by the components of a line into two parallel lines as 
equivalent line resistance. Ching [12] devised a method of transmission supplement charge allocation based 
on the topological analysis of power flows and DC linearized system model. This method is applied by 
assuming there is a line between the two ends of the shunt capacitor, which is exemplified as the fictitious 
load bus. Kaigue Xie [13] proposed an analytical model and algorithm based on the power flow, 
using Extended Incidence Matrix (EIM). In this method, firstly, a construction matrix EIM is defined as 
Matrix A, followed by defining the lossless power network (LPN) that provides the information about bus 
incidence relations and direction of active power in branch flow. Y.A. Altaturki [14] proposed the use of 
Current Adjustment Factor (CAF) method to allocate real and reactive power due to system linearity. In ref. 
[15] S. Nojeng, proposed determine of power flow based on the power factor aprroach. The result of this 
method had proven its validity and consistency. However, this method failed to consider the power factor of 
the load, though it is one of the factors associated with the transmission transfer capability. In general, 
the power factor does not change the active power, but it can affect power quality. 

The impact of the load characteristic (for instance; industries load, etc.) depends on the power flow 
contribution of each generator to the load in each transmission line. However, several methods have been 
used to the power flow tracing without the consideration of load characteristic. 

This paper is organized as follows: Section II describes the current Distribution factor method, while 
Section III describes the proposed method based on the load power factor. Application of the proposed 
method on the IEEE-30 bus system and discussion are presented in Section IV. Finally, Section VI concludes 
this study paper. 


2. DISTRIBUTION FACTOR METHOD 

The Distribution factor uses the A factor to determine the power flow in each line, through a DC 
line flow approximation. This principle uses the topology factor distribution by implementing a matrix of 
Generation Shift Distribution Factor (GSDF), factors Generalize Generation Distribution Factor (GGDF) and 
Generalized Load Distribution Factor(GLDF) to determine the contribution of the generator or the load on 
each line. 

A Factors are defined as a sensitivity factor, calculated as follows: 


— Gin Ze) 


1 () 


Aix 


Matrix A also is defined as Generalization Shift Distribution Factor (GSDF). 
Where: 
Zjx = Reactance of line from bus i to user k 
Zjx = Reactance of line from bus j to user k 
Xj; = Reactance of line from bus i to bus j 
For determining the load contribution to the line flow, the Generalized Load Distribution Factor 
(C factor) can be used, as follows: 


LF, a DiA-ik L, 
Cie) = x a = 
k 








(2) 


Equation (2) can also be defined as a C factor, where Fi-j is power flow in line i-j, Lk is load-k, R is 
Reference bus and k as bus-k. The Distribution Factor in the existing methods, expressed as Equation (1) 
to (2), only shows the contribution of the power flow to each generator and load based on the line impedance. 
This method, however, only considers the line parameters (line reactance), the power injected by the 
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generators and the load consumption, without considering the characteristics of each load. In reality, it does 
not reflect the characteristics of the load (ie. load power factor) which can affect the transmission line usage. 


3. PROPOSED METHOD 
Poor power factor will cause a decrease in transmission power transfer capability. Figure 1 shows a 
single line diagram, with the resistance ignored. 


P+AP 
x, 
Figure 1. Single line diagram 


The line impedance can be expressed as follows: 


Z,=R,+ JX, (3) 
The real power P can be expressed as follows: 


Where: 
= Angle betwen the reactive component with active componen in line where V and I are the voltage and 
@ current, respectively. 
Z, = Line impedance 

With assumption that X1 is line reactance of transmission line based on reference power factors, i.e 
(cos phi=0.8). Hence, by assuming RI=0, we can obtain: 


_ 172 
P=V~.X,.cos¢ (5) 


Assuming that the Vb (old) and Vb (new) are voltages in bus as constant before and after improving 
the load power factor, the following Equation is formulated: 


PiaX ota, = PiywX 


new~~ 1,new (6) 
Therefore: 
xX oon X ota 
I,new P 
new 
P 
old (7) 


With consideration of value from Figure 1, Equation 7 can be expressed as follows: 


X ota X oa 


l,new [ata | a. p22] 
Pia Pua (8) 


After eliminating the power flow in Piney) and Pyoig), Equation 8 can be written as: 


xX 
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By substituting Equation (9) to (1), we can obtain a new A factor as follows: 
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Meanwhile, the load contribution based on power factor for each line can be expressed as follows: 
Fa a DAs (new L, 
Chir (new = SL 
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To illustrate how the modified distribution factors can be used to determine the contribution of the 
generator and load to the line flow, a 3 bus system was considered for application, as shown in Figure 2. 


15.19 Mw MW. 1.65 MW 





13.54 Mw MW 


Figure 2. The 3 bus systems simple 


Case: This method was verified for the 3 bus samples, which consisted of one generator (G1) and 
two load (L2 and L3). GSDF and modified GSDF were obtained by using Equation (10 ) and (11), as shown 
in Table 1. The GSDF was then compared with the proposed GGDF method, based on load characteristics. 
These load characteristics were estimated based on the power factor condition in load 2 (by assuming GSDF 
having a power factor 0.8, taken as reference power factor). 

Similarly, to obtain modified GLDF, Equation (10) and (11) had been used to obtain the original 
GLDF and modified GLDF, as shown in Table 3 (by assuming GLDF to have power factor of 0.8, taken as 
reference power factor). 


Table 1. Comparison Between GSDF and Modified GSDF Based on the Power Factor in Load 2. 








Modified GSDF 
GSDF existing cos®=0.85 (lagging) cos®=0.90(lagging) 
0 -0.500 -0.250 0 -0.515 -0.274 0 -0.543 -0.319 
0 -0.500 -0.750 0 -0.484 -0.726 0 -0.456 -0.680 
00.500 -0.250 0 0.484 -0.274 00.456 -0.319 
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Table 2. Comparison Between GLDF and Modified GLDF Based on the Power Factor in Load 2 








ave Modified GLDF 
Grp Ressane cos ® =0.85 lagging cos ® =0.90 lagging 
00.500 0.250 00.515 0.274 0 0.543 0.319 
00.500 0.750 0 0.484 0.726 0 0.456 0.680 
0 -0.500 0.250 0 -0.484 0.274 0 -0.456 0.319 





Table 3. Contribution of Load in Each Line 
Contribution of load in each line (MW) 








Load Line GLDF Modified GLDF 
From To Method existing cos ® =0.85 lagging cos ® =0.90 lagging 
122 0.825 0.845 0.896 
L2 3-2 0.825 0.805 0.754 
Total 1.650 1.650 1.650 
1-3 10.155 9.830 9.214 
L3 2-3 3.709 3.709 4.325 
Total 13.540 13.540 13.540 





4. CASE STUDY 

In this case study, the proposed method was tested on the IEEE 30-bus system, as shown in 
Figure 3. The transmission network data used in the IEEE 30-bus system was referred from [14]. In this 
system, six generators were employed at buses 1, 2, 13, 22, 26 and 27, while 20 loads were located at buses 
2,3, 4, 7, 8, 10, 12, 14, 15, 16, 17, 18, 19, 20, 21, 23, 24, 26, 29 and 30. The generators served a total system 
load of 259 MW. 
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Figure 3. IEEE 30-bus system 


Table 4 shows the contribution of each load in line, obtained by using modified Distribution Factor 
method. As can be seen from the results, both methods showed difference in the power flow in line. 
The GLDF modified showed that the load contribution in line was in accordance with the changes of the 
power factor at load LS. For instance, by using the original GLDF method, load L5, contributed 46.673 MW 
to line 1-2. On the other hand, by using the proposed method, the contributions of load L5 were 50.016 MW 
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(cos ® =0.85 lagging) and 55.462 (cos ® =0.90 lagging), respectively. As load L7 was related to bus 5, there 
were also changes due to the connected line 5-7. Similarly, load L7 contributed 6.6675 MW (cos ® =0.85 
lagging) and 7.4923 MW (cos ® =0.90 lagging) to line 5-7. Both methods had the same total flow 
contribution of the lines related to the respective load. For example, the total flow contributed by L5 to line 
2-5, and 7-5 for both methods was 94.2 MW. This happened similarly to L7 and L30, which produced 22.8 
MW and 10.6 MW, respectively. 


Table 4. Contribution of generator Using the GLDF and Modified GLDF Based on 
the Power Factor in Load 5 





Contribution of loads L5,L7 and L30 in each line (MW) 





Load Line GLDF Modified GLDF 
From- To ee cos ® =0.85 lagging cos ® =0.90 lagging 

LS 2-5 48.67 50.02 55.46 
7-5 45.53 44.18 38.74 
Total 94.20 94.20 94.20 
L7 5-7 6.07 6.67 749 
6-7 16.73 16.13 15.31 
Total 22.8 22.80 22.80 
L30 30 - 27 6.24 6.24 6.24 
30 - 29 4.33 4.33 4.33 
Total 10.67 10.67 10.67 





5. CONCLUSION 

This paper has presented the ability of the modified Distribution Factors method for identifying the 
contribution of generator and load to transmission line flow, due to various load power factors. 
Determination of load contribution to electricity grid system has become important to the deregulated power 
system, especially for bilateral transaction in the electrical energy supply. By considering power factor 
variation from the users, the most efficient policy for determining the transmission charges for users (for 
instance: industry consumer) can be established. This case study results, based on IEEE 30-bus system, 
have shown that the proposed method can effectively determine the generator and load share in a fair manner. 
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